Intensity Interferometry with Anyons 
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A pairwise correlation function in relative momentum space is discussed as a tool to characterize 
the properties of an incoherent source of non-interacting Abelian anyons. This is analogous to the 
Hanbury-Brown Twiss effect for particles with fractional statistics in two dimensions. In particular, 
using a flux tube model for anyons, the effects of source shape and quantum statistics on a two- 
particle correlation function are examined. Such a tool may prove useful in the context of quantum 
computing and other experimental applications where studying anyon sources are of interest. 
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Anyons are particles that exhibit fractional quantum 
statistics in two dimensions. Existing somewhere be- 
tween bosons and fermions, anyons have been applied 
theoretically to a variety of problems such as the frac- 
tional quantum Hall effect (FQHE) I", "5| , high temper- 
ature superconductivity 3] , supersynimetry ^ |^ , and 
fault-tolerant quantum computing Experimentally, 
the quasi-particle excitations seen in the FQHE have 
been shown to possess anyon properties 0, ■ 

Fractional spin in two dimensions directly addresses 
the topological interpretation of quantum statistics. For 
identical particles in three or more dimensions, the a_ 
propriate homotopy group is the permutation group 
and there are only two homotopy classes. These corre- 
spond directly to fermions and bosons. 

However, in two dimensions ,^he corresponding homo- 
topy group is the braid group Q and an infinite number 
of homotopy classes are possible. As particles are ex- 
changed in relative position, each successive winding of 
the particles around each other cannot be smoothly de- 
formed into a finite number of configurations as is the 
case in higher (>2) dimensions. Each winding corre- 
sponds to a separate homotopy class which in turn can be 
categorized as a different classification of quantum statis- 
tics. The term anyon was coined by Wilczek 0, 0| 
to describe objects having such topological properties, a 
physical example of which is a charged particle bound to 
a magnetic flux tube in two dimensions. 

This relationship between topology and statistics was 
glimpsed by Aharonov and Bohm 'l2i and later realized, 
formalized, and expanded by others 10, 11, 13, 14, 15, 16, 
Il7| . Because the homotopy group is the braid group, this 
makes many multi-particle (>2) anyon problems nearly 
intractable. An essential analytical difficulty arises be- 
cause, unlike the case with the permutation group, the 
multi-particle wave functions for anyons cannot in gen- 
eral be written in a simple way in terms of the single- 
particle wave functions. Multi-particle anyon wave func- 
tions are, in effect, permanently entangled. However, it 
is exactly this topological property of anyons that makes 
them robust against decoherence and thus desirable can- 
didates for qubits in fault-tolerant quantum computing 



6]. 

Because quantum interference effects can be sensitive 
to quantum statistics, it is natural to inquire about 
the role of anyons in this context. As mentioned, the 
Aharonov-Bohm effect was an early probe into this fasci- 
nating problem. More recently, first order interference 
effects in a Mach-Zender-style interferometer for non- 
Abelian anyons in the presence of Aharonov-Bohm flux 
sources have been derived from applications of the braid 
group [l3|. 

In this work, intensity interferometry in momentum 
space, a second-order interference effect, is explored for 
incoherently emitted, non-interacting Abelian anyons. 
This is conceptually related to the recent q-Bose gas in- 
terferometry approaches for pions in heavy ion collisions 
[l^ |20|. However, q-bosons and q- fermions are not re- 
stricted to two dimensions. In that case, the generaliza- 
tion of quantum statistics is achieved somewhat differ- 
ently than for anyons |2ll |. 

Intensity interferometry is a useful tool for characteriz- 
ing the source geometry of incoherently emitted particles. 
In addition, the interference effect, expressed in terms of 
a second-order correlation function, might be regarded as 
an entanglement measure. Such a tool may prove useful 
in the context of quantum computing and other experi- 
mental applications where studying anyon sources are of 
interest. 

Intensity interferometry was originally developed by 
Hanbury-Brown and Twiss (HBT) as an alternative to 
Michelson interferometry to measure the angular sizes of 
stars in radio astronomy .22] . By correlating intensities 
rather than amplitudes, the measurement is insensitive 
to high frequency fluctuations that would normally make 
Michelson interferometry prohibitive. When applied to 
classical waves, the HBT effect is essentially a beat phe- 
nomenon. 

The field of modern quantum optics was spawned when 
intensity interferometry was quantum mechanically ap- 
plied to photons rather than classical waves. Now, two- 
and multi-photon effects are routinely studied in photon- 
ics. Two-fermion HBT in 2D condensed matter systems 
has also been reported |2^. The HBT effect was inde- 
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pendently applied to pions and other particles in elemen- 
tary particle physics in momentum space and is often 
called the Goldhaber-Goldhaber-Lee-Pais (GGLP) effect 
in that context [illilE^. 

Quantum mechanically, the two-particle correlation 
function is a measure of the degree of independence be- 
tween a joint measurement of two particles (1 and 2 be- 
low). Formally, the correlation function can be written: 



C. = 



Tr[patataa] -P(l,2) 
Tr[pata]Tr[pata] ^ P(1)P(2) 



(1) 



where p is the density matrix and and a are the cre- 
ation and annihilation operators for the quanta associ- 
ated with the appropriate fields of interest. When nor- 
malized to the single particle distributions as shown, C2 
is proportional to the relative probability for a joint two- 
particle measurement as compared to two single-particle 
measurements. If the measurements are independent, 
then C2=l. If the measurements are correlated, C2 de- 
viates from unity. In principle, C2 can be expressed as a 
joint measurement as a function of any degree of freedom 
of interest. In this treatment, we will focus on momen- 
tum space correlations. 

The second-order correlation function in momentum 
space is a powerful tool to probe several important prop- 
erties of a system. C2 is sensitive to the quantum statis- 
tics of the particles as expressed in the commutation re- 
lations for the field operators in Eq. (Q). Also, as de- 
termined by the form of the density matrix, C2 contains 
information about the space-time distribution of the par- 
ticle source in phase space as well as the pairwise inter- 
action and the quantum field configuration. 

Non-interacting identical particles can exhibit strong 
correlations, no correlations, or even anti-correlations de- 
pending on the specific field configuration. Possible field 
configurations for identical bosons include thermal states, 
coherent states, Fock states, and squeezed states H^]- For 
example, an incoherent thermal field configuration gives 
C2 = 2 for spinless bosons of equal momentum. In con- 
trast, a coherent source of bosons, such as laser light well 
above the lasing threshold, gives a constant C2 = 1 for 
all relative momenta. 

Geometric information about the source is contained in 
the shape of the correlation function as a function of the 
pair's relative momentum. The correlation function typ- 
ically approaches unity for large momentum differences. 
For incoherent sources, the relationship between source 
width and correlation width is essentially limited by the 
Heisenberg uncertainty principle such that 6q ~ 1/R 
where 8q is the width of C2 in relative momentum space 
and 1/i? is some characteristic width of the incoherent 
source generating the pairs. 

Keeping in mind the above discussion, this paper ex- 
amines the behavior of C2 in momentum space for a non- 
expanding, non-relativistic, incoherent source of anyons. 



An approximation to Eq. is used to highlight some 
features of intensity interferometry in momentum space. 
The quantum statistics are "tuned" , and, given some sim- 
ple source functions, the corresponding correlation func- 
tion is extracted. 

In momentum space for an incoherent source, Eg. ([T]) 
reduces to the Koonin-Pratt equation |23, 123, Oil 

C2(q)-1 = j dvK{q,r)S{v) (2) 
where the integration kernel is given by 



X(q,r) = |$q(r)|2-1. 



(3) 



The function $q(r) is the two-particle wave function 
in the center of mass frame of the pair, where q — 
(Pi ~ P2)/2 is the relative momentum and r = ri — r2 is 
the relative separation in that frame. The center of mass 
motion of the pair is not considered in this treatment. 
The source function, iS'(r), is the normalized probabil- 
ity distribution of emitting a particle pair with relative 
separation r. The integral is over the entire relative sep- 
aration space. In this context, incoherent means that the 
particles are emitted from the source randomly and in- 
dependently. Moreover, any potential time dependence 
of the source function in Eq. (0) has been integrated out 
as in Ref . . This integration eliminates the ability to 
recover emission-time ordering information between the 
pairs and places the focus on spatial information only. 

Using a flux tube model for anyons, a dynamical ap- 
proach is used to obtain two-particle non-relativistic wave 
functions for anyons in two dimensions 
this model, the quantum statistics are enforced through 
an interaction, the strength of which enforces the effective 
quantum statistics. In this picture, each anyon is anal- 
ogous to a charged particle bound to a tightly bundled 
magnetic field in two-dimensions. The analogy is not 
strictly one-to-one in this treatment because the scalar 
Coulomb potential between the charged pairs is not con- 
sidered; indeed, the system need not be electromagnetic. 
Nevertheless, the analogy paints a helpful physical pic- 
ture. 

The tightly bundled magnetic field associated with 
each anyon, B = V x A, points perpendicular to the plane 
of motion. We can select a gauge such that the vector po- 
tential, A, has the form Ar = and A^ = ^/27rr, where ^ 
is the magnetic flux through the plane. The particles are 
only permitted in regions of space where the magnetic 
field of the other particle is zero, so the particles exert no 
force on each other. Quantum mechanically, the interac- 
tion term takes the form of a minimal coupling such that 
the momentum operator takes the form p — > (— i V — eA) . 
While under the conditions described above, this "inter- 
action" exerts no force, it does adjust the phase of the 
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wave functions, permitting interference effects. The rela- 
tive wave function of the pair can then be used in Eq. jSJ, 
in combination with a normahzed source distribution, to 
obtain C2. 

In the center-of-mass frame, using the Hamihonian 



H 



Pi , {P4>-^f 



(4) 



the wave function is obtained for the relative motion of 
two free anyons in two-dimensional polar coordinates. In 
Eq. m is the mass of one particle (mi = TO2 — m) 
and a represents the anyon parameter with Q < a < \, 
where a = 0(1) corresponds to bosons (fermions). In 
this standard treatment, all of the details of the vector 
potential, magnetic flux, and gauge choice are included 
in the tunable parameter a. In Eq. the reduced mass 
of the identical pair, ^ = 17111712 / {mi + 7712) = m/2, has 
already been substituted. In the two-dimensional relative 
coordinate system of the identical pairs, the x-axis has 
been chosen along the direction of q. 

Applying Eq. Q to the time-independent Schrodinger 
equation, Hip = Eip, gives 



[( 



r dr ' 



1 d 



+ + q^]ijj = (5) 



where — mE is the square of the relative momentum of 
the pair (i.e. the energy is that of a free relative particle) 
and (t> is the angle between q and r. Eq. (0) admits 
solutions of the form 



V2 



E 



(^)^'''^^J\l~o.\{qr)e^"^ (9) 



00 (even) 



where / is the orbital angular momentum quantum num- 
ber. The sum is over partial waves of Bessel functions 
of the first kind of order |^ — aj where I is necessarily an 
integer, but a is not. The subscript q has been added to 
the wave function at this stage to remind the reader that 
we will be integrating over r and (j) using Eqs. (O and 
(O to obtain a function of q. 

Treating the initial pair wave function in Eq. (O as 
that of spinless bosons (the "boson basis" ) simplifies the 
problem of obtaining the appropriate exchange symmetry 
for the final two-particle wave function shown in Eq. ||^. 
In constructing Eq. lO, the partial wave expansion using 
Eq. lO is symmetrized under the exchange — > (0 — 
tt). Tuning the anyon parameter will result in a wave 
function that possesses the correct exchange symmetry 
for identical indistinguishable pairs. For the fermion and 
boson cases, Eq. (jSJ, with the help of Eq. JHJ, reduces to 
Eq. 0. But for < a < 1, Eq. © is the non-trivial 
two-anyon free particle wave function. 

Equation ||2J) is used to obtain an expression of C2 for 
anyons. While intensity interferometry can be used to 
image three dimensional sources, for simplicity we will 
only consider angle-averaged sources that are a function 
of r rather than r. The integration kernel in Eq. Q 
is angle-averaged in two dimensions and then denoted 
Kq^ {q, r) such that 



V'(r,0) =7VJu_„|(gr)e'' 



(6) 



where iV is a normalization constant. We want solutions 
that correspond to the appropriate symmetric (antisym- 
metric) boson (fermion) free particle wave function in 
relative coordinates in the limit of a — 0(1). Specifically, 



'^a.=oAr,(l^) = -^{e'^-''±e-'^-n (7) 

where the a = case corresponds to the symmetric so- 
lution (+) while the a — I case is the antisymmetric 
solution (-). 

Keeping in mind the identity 



i"J„(gr)e 



in(p 



(8) 



n— — 00 



one can construct partial wave superpositions of Eq. © 
that give Eq. (71) for the limiting cases of a = 0,1. 
Namely, 



KoAq,r) = [— d0|$,^„(r,0)|^] - 1. (10) 

Inserting the expression for the relative wave function 
from Eq. © into Eq. H10|l and performing the angle- 
averaging in two dimensions gives 



^— +00 

KoAq,r) = [2 ^ | J|,_„|(gr)|2] - 1. (11) 

/ — — oo(cvcn) 

In the limiting case when a = (0, 1), Eq. (|11|1 simplifies 
to, with the help of Bessel function properties. 



Ko^^,Aq,r) = {-irJo{2qr) 



(12) 



which is, by construction, the same angle-averaged ker- 
nels one gets by using Eq. iQ directly. For the more 
general case when < a < 1, the sum in Eq. 1)11(1 is eval- 
uated numerically to an appropriately converging order. 

Two sample source functions are used to illustrate the 
effects of source shape and size on C2: 
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FIG. 1: versus gro from a normalized Gaussian source 

given by Eq. \V6l in two dimensions. The curves represent 
different values of the anyon parameter a. The top curve is 
the boson case (a = 0) and the bottom curve is the fermion 
case (a = 1). From top to bottom, a =Q, 0.1, 0.2, 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8, 0.9, and 1. The source function, plotted in 
arbitrary units, is shown inset versus r/ro. 



Ss{r) = 



e(ro - r) 



(13) 



(14) 



They represent incoherent pair emission distributions 
normalized in two dimensions, 



POO 

2tt / rdrS{r) — 1, 
Jo 



(15) 



and where the effective source width is given by rp. The 
forms of the functions were chosen as examples to repre- 
sent typical localized incoherent sources. 

Integrating Eq. |(5J over the azimuthal angle in two 
dimensions using Eq. pi|l . and choosing either Eqs. Ijl^fl 
or (|14|l as the source function, gives an angle-averaged 
expression for C2 parameterized by a: 



C2„ (q) - 1 - 27r / rdrKo^ (q, r)S{r). 



(16) 



Figure^shows for an incoherent Gaussian source, 
Eq. lfO|l . Various values of a are displayed. The fig- 
ures are shown plotted versus the quantity gro, which 
accounts for the shifting momentum scale upon changing 
the source size. For example, as the width of the source 
decreases, the correlation function widens and vice versa. 

For the boson and fermion cases, the standard HBT re- 
sults for incoherent sources are recovered. The top curve 
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FIG. 2: from a normalized step function source given 

by Eq. 11411 . The anyon parameterization is the same as in 
Fig-Q The source function, plotted in arbitrary units, is inset 
plotted versus r/ro. 



in Fig. n represents the spin zero case and approaches a 
value of 6*2^^0 = 2 as g ^ 0. Here, it is more likely to 
measure a pair of non-interacting identical bosons in a 
state of zero relative momentum than to independently 
measure each particle with their respective momentum. 
That is, in the joint measurement, the bosons are not 
independent. This effect is a function of q and is the 
result of the Bose-Einstein statistics and the symmetric 
nature of the wave function. Similarly, the lowest curve 
represents non-interacting indistinguishable fermions and 
displays an anti-correlation in the joint measurement at 
low relative momentum, where C2„^i = as g ^ 0. Sim- 
ilar to the boson case, this is also only due to the quan- 
tum statistics and indicates again that a measurement 
of one particle affects the other quantum mechanically. 
As the relative momentum increases, the value of C2„ 
approaches unity, indicating that the joint measurement 
becomes less correlated. 

Starting from the top boson curve and moving down- 
ward are the correlation curves for various values of the 
anyon parameter. Like fermions, the anyons display an 
exclusion principle and always anti-correlate in the limit 
as 5 — > 0. However, other than the exclusion at g = 0, 
there is a natural trend of decreasing correlation as one 
interpolates between bosons and fermions for values of 
(7 > 0. The correlation function also approaches unity 
for large values of q. 

An example of the sensitivity of C2„ to the source 
shape is shown in Fig. El The box profile with a sharp 
edge introduces a ringing structure into the correlation 
function. 

Experimentally, the measured correlation function pro- 
vides information about the size and shape of the emis- 
sion source. Measured correlation functions can be in- 
verted to determine actual source distributions using 
methods like those discussed in Ref. (s^l- Also, informa- 
tion about anyon pairwise interactions can be extracted 
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by comparing measured correlation functions against pre- 
dictions made using hypothesized interactions in the 
Hamiltonian. Moreover, the degree of chaoticity of the 
source, testing the independent emission assumption, can 
be extracted by looking at the behavior of the correlation 
at low relative momentum. Finally, the correlation func- 
tion itself provides information about the degree of entan- 
glement between pairs versus their relative momentum, a 
potentially useful tool in the field of quantum computing. 

In summary, intensity interferometry in momentum 
space provides a tool that can be used to experimen- 
tally study the properties of anyons and their sources 
in two dimensions. Further theoretical studies using 
this framework can be pursued exploring various proper- 
ties of multi-anyon systems including anyon interactions, 
non-Abelian anyon correlations, and higher-order multi- 
anyon correlations. 
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